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Abstract 

Within the microscopic model based on the algebraic version of the resonating group method the role of the Pauli 
principle in the formation of continuum wave function of nuclear systems composed of three identical s-clusters has 
been investigated. Emphasis is placed upon the study of the exchange effects contained in the genuine three-cluster 
norm kernel. Three- fermion, three-boson, three-dineutron {'id') and 3q systems are considered in detail. Simple 
analytical method of constructing the norm kernel for Set system is suggested. The Pauli-allowed basis functions 
for the 3a and 3d' systems are given in an explicit form and asymptotic behavior of these functions is established. 
Complete classification of the eigenfunctions and the eigenvalues of the ^^C norm kernel by the *Be= a+a eigenvalues 
has been given for the first time. Spectrum of the ^^C norm kernel is compared to that of the system. 
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1. Introduction 

The question of the role of the Pauh principle in three-cluster systems goes back a long way in history. The 
^^G nucleus has been established to exhibit threc-alpha-cluster structure forty years ago. Since then many 
microscopic, macroscopic and semi-microscopic cluster models have been applied to analyze the structure 
of the ground and excited states of this nucleus. In particular, microscopic 3a calculation was performed 
within the resonating group method (RGM) by Kamimura iL| and within the generator coordinate method 
(GCM) by Uegaki [5]. Both calculations give reasonable results for the ground state of the ^^C and some 
excited states. However, wave functions provided by these models are very complicated and heavy to handle. 
Furthermore, although RGM ensures correct account of nucleon exchange between different clusters, an 
antisymmetry requirement on the total wave function can be violated by the improper truncation of model 
space. The measure of this violation is lacking. For example, in Ref. [T] the RGM calculation was performed 
with truncation to a space where the angular momentum I oi a — a relative wave function was fixed to be 
equal to zero. At the same time, even the lowest Pauli-allowed state of the ^^C represents the mixture of 
Z = 0, 2 and 4. This raises the question as to whether such truncation is consistent with the requirements of 
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the Pauli principle. As for the GCM, seven-dimensional numerical integrals are involved in the calculation 
and the integration over generator coordinate is replaced into summation over some mesh-points. However, 
the generator coordinate is chosen to be real and its domain is not well-defined, while only complex generator 
parameters ensure the existence of inverse transition from the generator parameter space to the coordinate 
space [3]. 

A number of macroscopic models were also applied to studying the ^^C nucleus (see, for example, Refs. 
[4l5p . Within such models, clusters are considered to be structureless particles interacting via local potentials, 
which reproduce experimental a — a phase shifts. But the antisymmetrization prevents a — a potential from 
being local. Furthermore, an accurate elimination of the Pauli-forbidden states is needed to get any reliable 
results on the structure of low- lying states of ^^C nucleus. By this reason, all macroscopic models fail to 
reproduce main properties of the 3a system. 

Orthogonality condition model (OCM) has been proposed in Ref. ^ as an approximation of resonating 
group method and applied to the investigation of the three-alpha structure in ■^'^C in Refs. [718] . In OCM 
the inter-cluster wave function is required to be orthogonal to the forbidden states and, consequently, the 
latter states are completely removed from the wave function. In two-cluster systems, OCM allows one to 
obtain only eigenfunctions of the antisymmetrization operator, not eigenvalues, which enter the Schrodinger 
equation and thus affect the dynamics of cluster system. Hence, in OCM a part of exchange effects is missing 
while it could be essential, especially in three-cluster systems. This is supported by the fact that the lack 
of binding energy of ^^C is observed in OCM compared to RGM calculations [9 . In Ref. [10 Horiuchi 
assumed that the Sa redundant solutions are all reducible to a — a redundancy and the Pauli-allowed 3a 
wave functions can be defined from the requirement of their orthogonality to the a — a forbidden states. 
Horiuchi also noticed that the Pauli-allowed basis functions should be classified by the number of total 
oscillator quanta and indices of the SU(3) symmetry (A,/i), because the second-order Casimir operator of 
the SU(3) group commutes with the operator of permutation of the nucleon position vectors. However, most 
of the three-cluster states are SU(3)-degenerate and the degree of SU(3)-degeneration increases drastically 
with increasing the number of oscillator quanta. Pauli-allowed basis functions corresponding to the same 
SU(3)-symmetry indices differs in their eigenvalues and should be labeled by additional quantum number. 
But construction method of OCM does not specify the character of the latter quantum number and, hence, 
in the case of SU(3)-degeneration neither eigenvalues nor eigenfunctions of the antisymmetrization operator 
can be uniquely determined. 

Recently Fujiwara and co-authors proposed a new type of three-cluster equation that employs a two- 
cluster RGM kernel for the inter-cluster interaction |!lF. The authors applied this equation to systems 
composed of three dineutrons and three alpha-clusters. Results for binding energy for 3a system, obtained 
by diagonalization of the Hamiltonian using [3] symmetric translationary-invariant harmonic oscillator basis, 
have been compared to the RGM and OCM calculations. Energy of the ^^C ground state given in Ref. [11] is 
closer to experimental data than the value provided by OCM, but still underbound by 1.5 MeV compared to 
the 3a RGM. The reason is that the three-cluster equation used in Ref. ||1T] does not include some exchange 
effects contained in the genuine three-cluster norm kernel. 

The three-cluster norm kernel contains complete information about the Pauli-allowed model space for the 
relative motion of clusters. Eigenvalues and eigenfunctions of the norm kernel are solely determined by the 
assumed internal wave functions of the clusters. Hence, careful analysis of these quantities is very useful for 
understanding the structure of three-cluster systems. However, there are few papers concerning this question. 
Eigenvalues and eigenfunctions of the norm kernels for 3a and ^^0-t-2a systems have been calculated in Ref. 
[12) . They have been characterized by the SU(3)-indices and an additional quantum number, but the meaning 
of the latter number has not been established and no classification with respect to this number has been 
provided. The eigenstates of the norm kernel were expanded by the harmonic oscillator wave functions in 
the rectangular coordinate representation. Only three simplest Pauli-allowed states of the ^'^Mg have been 
given in Ref. [12 , while eigenfunction of the norm kernel for ^^C have not been presented. Coincidence 
of some eigenvalues of the ^^Mg with the eigenvalues of the two-body norm kernels of ^Be= a + a and 
20Ne= a-|-^^0 at large number of quanta has been noticed, but this property has not been taken advantage 
of in classification of the three-cluster Pauli-allowed states. 

Summarizing, there are several problems still outstanding; construction of eigenfunctions of the antisym- 
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metrization operator for the three-cluster systems in simple, tractable form, complete classification of these 
states (i.e., resolving SU(3)-degeneration problem) and elaboration of the truncation method of the model 
space consistent with the requirements of the Pauli exclusion principle. 

For the case of three-cluster systems composed of an s-cluster and two neutrons we have shown that all 
these problems can be resolved within the discrete representation of a complete basis of allowed states of 
the multiparticle harmonic oscillator (classified with the use of the SU(3) symmetry indices and defined in 
the Fock-Bargmann space) [13] ■ Classification of the eigenvalues of the three-cluster systems with the help 
of eigenvalues of the two-body subsystem was suggested in Ref . [T3' . We observed that asymptotic behavior 
of basis functions consistent with the requirements of the Pauli principle gives an indication of possible 
decay channels of a three-cluster nucleus and allows us to specify the most important decay channels. Such 
asymptotic behavior corresponds rather to the subsequent decay of the three-cluster system than to the 
so-called "democratic decay" associated with the hyperspherical harmonics, which are widely used for the 
description of three-cluster systems. 

In the present paper we extend our approach to the three-cluster systems composed of identical s-clusters. 
In Section[2]the theoretical model is briefiy explained. Norm kernels of three-fermion and three-boson systems 
are considered in detail in Section [3] and Section IH respectively. Section [5] is devoted to the analysis of the 
eigenvalues and eigenfunctions of the norm kernel of the threc-dineutron system. The Pauli-allowed states 
of the Sa-system are constructed and classified in Section [51 The asymptotic behavior of these functions 
is established and compared to the asymptotic form of the eigenfunctions of the "^H-f-n + n norm kernel. 
Concluding remarks are given in Section [71 The norm kernel of the 3a system defined in the Fock-Bargmann 
space is given in an explicit form in Appendix lAl 



2. Theoretical model 



The RGM wave function of three-cluster system is represented as an antisymmetrized product of internal 
cluster functions and wave function of the relative motion of clusters. The latter is found by solving the 
integro-differential equation, which is obtained from the Schrodinger equation upon integrating over single- 
particle coordinates. The internal cluster functions are considered to be fixed. Expanding the wave function 
of cluster relative motion into the basis of the Pauli-allowed multi-particle harmonic oscillator states \['„({r}), 
one can reduce the integro-differential equation to the set of linear algebraic equations for the expansion 
coefficients. Further simplification can be achieved by transition from the coordinate space to the space of 
complex generator parameters (Fock-Bargmann space [E]), because the Fock-Bargmann image \1/„({R}) of 
the Pauli-allowed state 4'„({r}) is much simpler than the original. Here, {r} is a set of nucleon position 
vectors, {R} is a set of generator parameters with the help of which the dynamics of the considered degrees 
of freedom will be reproduced, and n is a set of quantum numbers of basis functions. 

The complete basis of the Pauli-allowed harmonic oscillator states defined both in the coordinate space 
and in the Fock-Bargmann space is generated by the Slater determinant <I>({R}, {r}) composed of the single- 
particle orbitals. It is advantageous to use the modified Bloch-Brink orbitals as spatial wave functions of the 
nucleons, because each orbital serves as the kernel of the transformation from the coordinate representation 
to the Fock-Bargmann representation and, at the same time, it is generating function of the basis of a 
harmonic oscillator. For three-cluster systems composed of three s-clusters containing A/j < 4 nucleons, only 
orbitals of the following form are required: 



0fc(ri 



-ij exp (^-ir2 + y2(Rfc • r,) - ^R^^ , i & 



where complex generator parameters Rfc are independent variables of the wave function 'I'ndRfe}), which 
describe the positions of the clusters in the Fock-Bargmann space: 

^fc and ffk are vectors of coordinate and momentum, respectively. 
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A complete basis of Pauli-allowed states in the Fock-Bargmann space and their eigenvalues can be obtained 
from the overlap integral I{Sk, R/c) of the two Slater determinants, which are usually defined as follows: 



/(Sfc,Rfe) = J $(Sfc,{r})$(Rfc,{r})dT. 



Here, integration is performed over all single-particle variables. The motion of the center of mass is eliminated 
by transition from the vectors Rfc to the Jacobi vectors 

Re™ = ^ (AiRi + A2R2 + A3R3) , 



,jMA^f^^_A,^^±^\ b../^(R.^R3). 
V A V A2 + A3 J ' V ^2 + A3 ' 

Then 

/(Sfe, Rfc) = exp(Rcm ■ Scm) /(a, b; a, b), 

where /(a, b; a, b) is a translation-invariant overlap integral generally known as the norm kernel. 
The Pauli-allowed basis states vE'„(a, b) are the eigenf unctions of the norm kernel: 

A„^'„(a,b)= J /(a,b;a,b)*„(a*,b*)d/i~ ^, 
while A„ are its eigenvalues. The Bargmann measure dfiii \3 is defined as follows: 

d^i^,b = exp{-(aa )} exp{-(bb )} 

The Hilbert-Schmidt expansion of the norm kernel 
/(a, b; a, b) = ^ A„^'„(a, b)1'„(a, b) 

n 

can be interpreted also as density matrix of mixed system [14] , with A„ being the elements of the diagonalized 
density matrix, which are proportional to the realization probability of the system states defined by the 
corresponding eigenfunctions. 

Examine next the explicit expressions for the norm kernels of those three-cluster systems which comprised 
of three identical s-clusters, such as — n+n + n, three-boson system, ®n =^ n+'^n+'^n and ^^C=a + a + a. 



3. Norm kernel of three-fermion system 



In this section we consider norm kernel of three identical fermions (for example, three nucleons with the 
same spin and isospin projections), which is the main building block for the norm kernels of three-cluster 
systems composed of three identical clusters. 

Instead of Jacobi vectors a and b let us introduce new complex vectors 

a + ib „ a. ih 
A= B = —. 

Similar transformation of creation and annihilation operator has been shown to simplify essentially con- 
struction of a complete bispherical harmonic oscillator basis in three-nucleon system |16j . Further we shall 
demonstrate that such transformation facilitates also construction of SU(3)-basis for three-fermion system 
and, more importantly, solving eigenvalue and eigenfunction problem for the three-cluster norm kernel. 
In terms of vectors A and B the norm kernel for three identical nucleons has the following form: 
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/„+„+„ (A, B; A, B) = ^ exp { (AA) + (BB) } + ^ exp {e^^ (AA) + e"^^ (BB)} + 

+ ^ exp {e-^*(AA) + e'^ (BB)} - i exp {(BA) + (AB)} - 

- ^ exp {e*^ (BA) + e^^^ (AB)} - ^ exp |e-*^ (BA) + e** (AB)} . 

The latter expression can be considered as the result of action of the antisymmctrization operator A on 
the first exponent, which contains a compete basis of both the Pauli- allowed and Pauli- forbidden harmonic 
oscillator states: 

/„+„+„ (A, B; A, B) = liexp{(AA) + (BB)}, 

where A is defined in the Fock-Bargmann space as 

i/(A,B)==/(A,B) + /(e'*A,e-*^B) + / (e"'^ A, e'^B^ - 

- /(B, A) - / (e'*B, e-*^ a) - / f e'^^B, e'^ a) . (1) 



First term corresponds to the identity permutation, second and third terms describe two cyclic permutations, 
while the last three items are associated with pair permutations and, hence, have the minus sign. So, 
permutations in the system of three identical nucleons correspond to the rotation by 120 degrees in the 
plane going through all nucleons. Evidently, such rotation relates to the transition from one Jacobi tree to 
another: 

Ai=e*^B; Bi = e"'^ A; A2 = e'^^B; B2 = e'^A, 
where 

. _ai.2+?bi_2 ai.2-«bi^2 
^^^^ ■ V2 ' = ^ V2 • 

Vectors ai^2, ^1,2 belong to the alternative Jacobi trees. 

Expanding the norm kernel ^A, B; A, B^ in powers of vectors A and B we obtain two sums: 

00 

/„+„+„ (a, B; a, b) = ^ |(AA)'"(BB)™ - (BA)™(AB)™} + 

rn— 1 ^ ''^ 

00 00 ^ 

+ W TT TTT ((AA)'"+3'=(BB)" + (AA)'"(BB)'"+3'=- 

'^ f^ 2m\{m + 3k)W ' ^ ' \ 1 \ > 



(BA)"+3'=(AB)" - (BA)™(AB)"+3'^} 



Either of two sums generates the states with definite number of oscillator quanta v: the one-fold sum 
corresponds to v = 2m, while the two-fold sum contains the states with v — 2m + 3fc. It should be noticed 
that only expressions, which are antisymmetric with respect to nucleon permutations, enter the above 
expansion of the norm kernel. 

These sums each are the superpositions of partial norm kernels i^^ possessing definite SU(3)-symmetry 
and labeled by the SU(3) indices (A,/2) (see Ref. [T7]): 



2(m!) 



^ m 

— {(AA)'"(BB)'" - (BA)"(AB)"} = ^ ^(2„._2M,M) (a, B; A, b) , 



^1=0 



„ ^ ((AA)"+3'=(BB)" + (AA)™(BB)™+-'^'^ - (BA)"+^'=(AB)" - (BA)"VAB)"+3'^'| 

2m!(m + 3fc)! L J 



■ i(2m+3k-2p,,p.) ; 
/i=0 



where 



X {(AA)'"-'^(BB)'"-'^2i^i i-m + fi,-m + JX; -2m + 2fl; z) - 

- (-l)'^(BA)'"-'' (AB)'"-''2Fi (-m + /i, -m + /x; -2m + 2/2; i)} ; 



X {{(AA)'"+3'=-'^(BB)'"-'^ + (AA)'"-'^(BB)™+3*^-'^} x 
X 2-Fi (— m + /i, — m — Sfc + /(x; —2m — Sfc + 2/i; 2;) — 
-(-l)*^ |(BA)'"+3'=-'^(AB)'"-'^ + (BA)'"-'^(AB)'"+3'=-'^| x 
X 2F1 (— m + Ji, —m — 3k + fi; —2m — S/c + 2/2; i)} ; 

Here (— a;)„ is the Pochhammer symbol, defined as follows: 
(a; — my. 

2-Fi (a, (3; 7; is the hypergeometric function with variables 
([AB][AB]) ^_ ([BA][AB]) 



z = 



(AA)(BB) (BA)(AB) 

The number of oscillator quanta i> relates the SU(3)-symmetry indices (A, fi) as = A + 2/[t. 
Partial norm kernels ai'e normalized to the dimensionality of (A, p.) SU(3)-representation: 



/ 



(A n A RW» - (A + l)(/^ + l)(A + M + 2) 
«(A,/i) (A, B; A, B)o!/XA,B = 5 • (2) 



It is remarkable that «(2m-4/i,2M) = 0- Hence, the latter representation contains only the Pauli-forbidden 
states. Evidently, the state corresponding to zero number of oscillator quanta also can not be realized, i.e., 

i{o,Q) = 0. 

Basis functions characterized by zero total orbital moinontum L = arc generated by the i{2m+3k-2p,,p,) 
representations, provided that both SU(3)-indices are even: 2m+3fc = 21/, jl = 2fi. Explicitly, these functions 
are expressible in terms of hypergeometric functions 2-^1(0;, /?; 7; Z), with the variable 



A2B2 

Namely, 

2-Fi (—1/ + n + fi, —n + /u; — 1/ + 2/t + i; , u — 3k = 2n; 
(AB) 2-Fi {-V + n + 1 + fj,, -n + 11] -V + 2fj, + Z) , - 3fc = 2n + 1, 

6 



X 



^2(v + 3fc - 2/i)!(iy - 3fc - 2^)! 



(2z/ + 1 - 4^)! / {2v - 2pL)\\{2v - 4/i - 1)!!(2/^ - 1)!! 



y (2^)!(2i/ + 1 - 2/^)! y (2i/ - 2^ + 1)!!(2j/ - 4/i)!!(2^)!! ' 

A set of quantum numbers of functions Xi^^^'^'^^^j^^ (A, B) includes the total number of oscillator quanta 2v, 
the SU(3) symmetry indices (A,/x), and additional quantum number fc for the SU(3)-degenerated stateJ^. 
Total orbital momentum L and its projection M are equal to zero and, hence, will be omitted in what 
follows. 

Functions Xf^^^^"^^^) (-A, B) comply with requirements of the Pauli exclusion principle. It is easy to verify 
this acting on these functions by antisymmetrization operator A defined in ([1]). 



1, 1 < fc < [^]; 
0, A: 0. 



Here [a;] denotes an integer part of the number x. 

Hence, at each value of the number of oscillator quanta v there are [^^^^] Pauli-allowed basis functions 
and one Pauli-forbidden state belonging to the same SU(3) representation {2i' — 4fi, 2/i). All the Pauli-allowed 
states have unit eigenvalues, while the Pauli-forbidden states correspond to zero eigenvalues. This implies 
that a norm kernel of three-nucleon system can be treated as the density matrix of the pure system. Indeed, 
in this case dynamics of all degrees of freedom is considered and no averaging is performed. 

In representation of A and B vectors an SU(3)-projected norm kernel i(2j/_4^^2/j) (A, B; A, B) has a diag- 
onal form: 

[^] 

Z(2.-4.,2,)(A,B;A,B)= ^ xl^l';X,{A,B)xl,^^^^^^^^^^ 

k=l 

Hence, solving the eigenvalue and eigenvectors problem for three-nucleon system becomes a trivial procedure. 
It is also easy to distinguish basis states belonging to the SU(3)-degenerated representation by the additional 
quantum number fc, which explicitly enters the expression for the basis functions. However, due to the fact 
that all eigenvalues of the Pauli-allowed states equal to unity, eigcnfunctions of the norm kernel are not 
uniquely determined. Any unitary transformation applied to the functions X(2i''!^4^^2/t) (A-i B) gives a new 
set of basis functions, which are also antisymmetric with respect to the nucleon permutation. For example, 
hyperspherical harmonics can be successfully applied to studying three-nucleon systems. 



4. Norm kernel of three-boson system 



In this section we consider norm kernel of three identical bosons /3{,(A, B; A, B), which generates basis 
functions needed for the description of three-cluster systems composed of three identical boson clusters, such 
as n n n and ^^C=a + a + a. 

In terms of vectors A and B the norm kernel for three identical bosons is very similar to that of three 
fermions: 

/3b(A, B; A, B) = I exp {(AA) + (BB)} + ^ exp {e^*(AA) + e"'* (BB)} + 
+ 1 exp {e-'*(AA) + e'*(BB)} + ^ exp {(BA) + (AB)} + 
+ i exp {e^* (BA) + e"'* (AB)} + ^ exp {e"^* (BA) + e'* (AB)} . 

^ An upper index of basis functions eorresponds to the minimal of the number of quanta along A and B 
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The only difference between the norm kernels of three fermions and three bosons is that the latter is sym- 
metric with respect to nucleon permutation, while the former is antisymmetric. Naturally, /3b(A, B; A, B) 
can be considered as the result of action of the symmetrization operator P on the first exponent: 

hb{A, B; A, B) = Ipexp {(AA) + (BB)} , 
where P is defined in the Fock-Bargmann space as 

P/(A, B) = /(A, B) + / (e^^ A, g-^^b) + / (e"^* A, g'^b) + 

+ /(B, A) + / (e'*B, e-'^ a) + / (e-'^B, e'^ a) . (3) 



Meaning of each term is the same as in Eq. ([T]). 

Expansion of the norm kernel /3b (^A, B; A, B^ in powers of vectors A and B gives the result: 

(a,B; A,b) = {(AA)'"(BB)'" + (BA)'"(AB)'»} + 

OQ 00 ^ 

+ E E 2mHm + 3kV {(AA)™+^'=(BB)™ + (AA)"(BB)™+-^'^ + 

m=Ofe=l 

+ (BA)"+3fe(AB)" + (BA)"(AB)'"+^''^} . 

Of course, only expressions, which are symmetric with respect to boson permutation, enter the above expan- 
sion of the norm kernel. Notice, that the above-mentioned symmetry results in essential reduction of the basis 
states, which can be realized in three-boson norm kernel. Indeed, the first exponent exp |(AA) -I- (BB)|, 
which generates both the Pauli-allowed and the Pauli- forbidden basis functions, contains six times more 
states than /3b ^A, B; A, B^ norm kernel. 

By analogy with three-fermion system, each sum corresponds to the fixed number of quanta 9 and can be 
presented as a sum of partial norm kernels belonging to the SU(3)-representations (A, fl) : 



2(m!) 



— {(AA)"(BB)'" + (BA)"(AB)"} - Y.hm-2i.,i.) (a,B; A,b) , 



p.=0 



((AA)"+3'=(BB)™ + (AA)™(BB)™+''^'^ + (BA)™+3'=(AB)" + (BA)"(AB)™+3'=]. = 

i-l-3fc)! L J 



2m!(m -I- 3fc) 



where 

(-l)^(-m)p(-m)p 



*(2m+3fe-2;t,/i)) 

M=0 



X |(AA)™-^'(BB)™-''2i^i {-m + fl,-m + fi; -2m -I- 2/1; z) + 
+ (-l)^'(BA)™-^'(AB)™-^'2Fi i-m + fl, -m + fl; -2m + 2/2; i)} 



*(2m+3fc-2p,/3) 



(a B a b]- (-iy''(-m)^(-m-3% (rABUABir 

^A, t!. A, a J - 2^,^^ ^ 3fc)!/i!(_2m - 3fc + /i - 1). ^^^^^ ^^^^^ 
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X ||(AA)™+''^'^-^'(BB)'"-'^ + (AA)"-'^(BB)"+^'=-^'| x 

X2^i {~rn + fl, — m — 3fc + /l; —2m — 3fc + 2/1; z) + 

+ (-1)'' |(BA)™+3'=-^'(AB)™-^ + (BA)'"-^'(AB)'"+3''^-^'| x 

X2F1 (— m + p., — m — 3k + ft; ~2m — 3fc + 2/1; z)} . 

Arguments z and z of the hypergeometric functions are the same as those in the previous section. 

In contrast to three-fermion system, i(2m-4/j,2/j) 7^ 0, while j(2m-4/i-2,2/i+i) — 0. Therefore, for three- 
boson system all the states belonging to {2iy — 4/i, 2fj.) representation are allowed by the Pauli principle. 
Basis functions characterized by zero total orbital momentum L = are of the following form: 



CA'2^)(A,B)=^(7/4X,) (A«.+B6.) A2„-2.B2„-2.[AB]2- X 

2F1 (— + n + II, —n + fi; —v + 2/t + i; , v — "ik ~ 2n; 
(AB) {-V + n + 1 + fi, -n + n; -ly + 2fi + ^; Z) , u - ik ^ 2t 



(4) 



^2{1 + dk,o){t^ + 3k- 2ii)\{v - 3fc - 2/i)! 



(2i/ + 1 - 4/i)! / {2v - 2ii)\\{2v - 4/i - l)!!(2/i - 1)!! 



y (2/i)!(2i/ + 1 - 2/i)! y (2z^-2/i+l)!!(2j/-4/i)!!(2/i)!! 
Functions '4'^(2y-i'^2 )(-A'-B) ^'^^ eigenfunctions of the symmetrization operator P defined in ([3|): 



V — 2^ 



Finally, at each value of the number of oscillator quanta [-^^-^j + 1 basis functions, belonging to SU(3)- 
representation (2z/ — 4/x, 2/i) and having zero orbital momentum, can be constructed, and all these states are 
consistent with the requirements of the Pauli principle. 

In the next section we will demonstrate that the Pauli-allowed states of three-cluster systems composed of 
identical boson clusters can be represented as a superposition of eigenfunctions of three- 

boson norm kernel. Indeed, the structure of the latter functions reflects the symmetry with respect to the 
permutation of clusters as a whole and, in addition, these functions are the eigenfunctions of the second-order 
Casimir operator. Further it will be shown that using ■0^^2^'^^~^2^j(A, B) functions we essentially simplify 
the eigenvalue and eigenfunction problem for three-cluster systems. 

5. Norm kernel of three dineutrons 



Norm kernel of three-dineutron system /3d' can be obtained by squaring of the three-neutron norm kernel 



hi, A,B;A,B -3 



/A B A B 



Replacement 
A 



-^A, B^^B, A-^^A, B^^B 
V2 %/2 V2 



V2 ' V2 ' V2" 

is required to ensure correct elimination of the center-of-mass motion in Sd' system. Such scaling trans- 
formation of the Jacobi vectors accounts for the fact that each cluster in the id! system consists of two 
neutrons. 
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The norm kernel I^d' can be presented as a three-termed antisymmetrization operator acting on the 
exponent, which contains a complete basis of both the Pauli-allowed and the Pauli-forbidden states: 



hd'iA, B; A, B) = + i2 + A3) cxp {(AA) + (BB)} 



Each of three above-mentioned operators combines certain group of permutations. 

The first operator Ai coincides with the symmetrization operator P of three-boson system defined in ([3]). 
Hence, Ai — P corresponds to the permutations of dineutron clusters as a whole. 

The second part of antisymmetrization operator is defined as 



^2/(A,B) = 2^/(^-^,-| 
B A 



+ .f[-e 



.A 



B 

"2 



/ 



+ / 



B 

Y' 



A 

y 



+ / 



B 

7' 



.A 

y 



(5) 



It is associated with permutations of nucleons between all three dineutron clusters and includes simultaneous 
pair or cyclic permutations of spin-up and spin-down neutrons. Notice that only the first term of operator 
A2 results in a rearrangement of clusters due to the permutation of identical neutrons. Remaining five terms 
can be obtained by permutation of clusters as a whole following the rearrangement associated with the first 
term. 

Finally, the last operator A3 is a little bit more complicated: 



A3^Al + Ai 



A 



3 ' 



where 



i^/(A,B) 



i+/(A,B) = -2 / 



-2|/( 



/A + B A + B 



2 ' 
A + B 



^, .2,A-|-B _2,A + B 
+ / e ' — 7, — ,e 3 — — 



A + B 



e'^A + e~'~B 



1 r 

2 



e'^A- 



-B 



e 3 



+ / 
+ / 

i3-/(A,B) = -2^/ 
+ / 
+ / 



-A 



^B 



-B 



e-— 



e»— A + e 
e"*^A + e'^B 
-'^A + e'^B 



e 3 



-A 



e'-A 



^B 



-B 



A + e'~B 
-*¥a + e'^B 



+ 



-A 



-B 



e 3 



-A 



-B 



(6) 



Operator A3 combines permutations of nucleons between two clusters, with the third cluster being a specta- 
tor. In this case the spin-up neutrons undergo cyclic permutation, while the spin-down neutrons experience 
pair permutation and vice versa. Depending on which cluster is a spectator, we obtain operators A3, or 
A'^ . The three-termed structure of each operator corresponds to three cyclic permutations of clusters as a 
whole. 

It is natural to seek the eigenfunctions of the id! norm kernel in the form of superposition of the eigenfunc- 
tions of the three-boson norm kernel, because the dineutron clusters should obey boson statistics. Hence, 
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the Pauli-allowed basis functions of the three-dineutron system ^{2v-4,ij.,2ij.)k{-^^ characterized by positive 
parity, zero values of total orbital momentum and spin L = S = and belonging to the {2p — 4/x, 
representation of the SU(3)-group can be written as: 

[^] 

*(2.-4.,2.). (A, B) = ^ V'(T.-T.%) (A, B). (7) 

g=0 

Expansion coefficients D'^2i^^-^'it^2tJ.)k' ^^o^S with eigenvalues ^{2v-4ij.,2ti)k corresponding to the Pauli-allowed 
states ^(2i/_4^^2^j)i. (A, B), can be found by diagonalization of the thrcc-dincutron system norm kernel 
*(2!^-4([i 2ij.) projected on the states with the fixed number of quanta and definite SU(3)-symmetry: 

i\i-A„,2>^)= J2 Yl (2;^,z^-3g-2/.|2!.,z.-35-3^-2/.)3d' X 

q=0 q=0 

^ \'f^{2,y-4n,2tj.)f^{2iy-4:fj.,2fi) '^^{2i^-4fj.,2fi) 'f^ {2u-4fi,2fj.) J ' 

In order to construct characteristic equation, let us consider the action of the antisymmetrization operator 
^ = Ai + ^2 + ^3 on the functions V'('2^!_^4^^2/i) • 
It is easy to verify that 

^^i^(T.X%)(A,B) = ^^^p4,!2.)(A,B) = ^V(2;X!^,)(A,B). 

Notice that all terms entering the expressions specifying operators Ai or A2 give the same result acting on 
basis functions V'^2"^!.%^2,«) • 

The last part of antisymmetrization operator provides nonzero result only for {2u, 0) representation: 



1 fi '''/^l"' 1 

--i3 <,"^V^o ^ (A, B) = 0— V VrV^VT^^^CA, B) , X 

(2^)! 



+ 3q + 3q)l{i^ -3q- 3g)!(:/ + 3g)!(z/ - 3g)! 
Consequently, the matrix elements of the secular equation look like 

(2., . - 3. - 2,12., . - 3. - 3,- - 2,)3. = (l + ^) ^.o - ^^4 ^(1 + ^(1 + .,-,) 



+ 3q + 3qy.{iy -3q- 3qy.{u + 3qy.{u - 3g)! " 



The action of the antisymmetrization operator A on the basis functions ^'(21/^4/2/^) results only in diagonal 
matrix elements of the characteristic equation for /x ^ 0: 



MT.-r.!2;)(A,B) = (1 + ^) ^(2;!.'^4;X)(A,B), , ^ 0. 



By this reason, the Pauli-allowed basis functions of the 3d' system belonging to the SU(3) representation 
{2v — 4/x, 2/Lt) with /x ^ coincide with the eigenfunctions of the three-boson norm kernel and have identical, 
although not equal to unity, eigenvalues A.(^2v-4ix,2ij.)k '■ 



*(2.-4^,2^). (A, B) = r^2~-4X) (A, B), A(2._4^,2^), = 1 + 0<k< 
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V — 2lJL 



The Pauli-allowed states of (2i/, 0) representations are characterized by the same eigenvalues at a given 
number of osciUator quanta v. However, there is also one Pauli-forbidden state, corresponding to zero 
eigenvalue: 



A(2.,o) - (1 - 4,0) {\ + , < fc < 



Eigenfunctions ^'(2,^,0)^ (A, B) of the antisymmetrization operator have more complicated form, being a 
superposition of functions "^^^^vvi) (^^^ ©)• However, with increasing the number of oscillator quanta 
the Pauli-forbidden state vE'(2^.o)fc=o(-^' takes simple analytical form: 



*(2.o)..o-^^— =(a +a, +a 



2 




«(A,B) 



It is quite obvious that such asymptotic behavior corresponds to the case when a dineutron moves away 
from the "^n —"^ n +^ n subsystem (tetraneutron), with the latter being in the Pauli-forbidden state. 

Similar behavior is peculiar also for the Pauli-allowed states. By means of the unitary transformation the 
eigenfunctions of the id! norm kernel corresponding to nonzero eigenvalues can be rearranged so that 

^'(2,.-V,2,.). ^ ^^2^-4^,2^). = ^ ('^(2r-t,2p) b) + (t\2u-t^.,2^.) bi) + 0(2^-^^,2^.) (^2, ba) 




jc\i\. -3q- 2^,)^l^l%%^ (A, B). (9) 



Here functions ^4 2 ) b) defined as 



^(2i/-4m,2^) b) = ^(2,y-4Ai.2M) I 



^ A^fjrf 2 1 [ab]2^a2''-2m-2M^j2m-2M 



1 \ [ahf 

X2F1 I —ly + m + fi, —TO + /i; -i^ + 2/i + - ; Z , Z = 



2' 7' a2b2 
with N'^^^I'i^ 2fj,) being a normalization coefficient: 



^2^-2^ ^ / (2Z. + 1 - 4/i)!(2z. - 2/i)!!(2t. - 4^. - l)!!(2/i - 1)!! 

(2^-4m,2m) y (2/i)!(2i/ 1 - 2/i)!(2i/ - 2to - 2/i)!(2m - 2^)!(2// - 2/i l)!!(2iy - 4/i)!!(2/i)!! ' 

— 1/2 

^2fe ^ 3(7 — 2/i) are the so-called Kravchuk polynomials of a discrete variable [TB] , which are specified and 
orthogonal on the interval 2^ < m < 2v ~ 2fj,: 

^1/2. 3 2/:)-PV'''' (2^-4/.)!2j^i(-2fc + 2M,-^ + 3g + 2/z;-2i. + 4/.;2) 
^''^ V2/ ^(2fc - 2/i)!(2i/ -2k- 2fi)\{iy - 3q - 2/i)!(j/ + 3g - 2/i)! ' ^ 

Hence, index of the Kravchuk polynomials serves as additional quantum number for the SU(3)-degenerated 
states. In Ref. [13, we have shown that a family of eigenfunctions of the antisymmetrization operator com- 
bining the states ^'(2i/-4;i,2/i)fc with the same index k corresponds to a certain binary decay channel of a 
three-cluster system into a two-cluster subsystem occurring in a ground or an excited harmonic-oscillator 
state and a remaining cluster. In turn, eigenvalues A(2,y_4p,2/j)fc of three-cluster system tend to the eigenval- 
ues of a two-cluster subsystem with increasing the number of oscillator quanta v. However, all eigenvalues of 
the Pauli-allowed states of "^n 4-^ n subsystem are equal to unity, in contrast to the eigenvalues of the '^H+n 
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subsystem of the '"'H considered in [13j . By this reason, eigenvalues A(2iy_4^,2/i)fc of the 3d' system are equal 
and tend to unity for 1 < A: < [^^-^] • Hence, strictly speaking, the Pauli-allowed states of the 3d' system are 
not uniquely determined and any unitary transformation can be applied to these states without violating the 
Pauli principle. It would seem that a basis of hyperspherical harmonics also could be used for study of the 
3d' system. However, existence of a Pauli-forbidden state makes this statement quite doubtful. Indeed, since 
a Pauli-forbidden state and the Pauli-allowed states have different eigenvalues, no unitary transformation 
can leave the diagonal form of the norm kernel intact. Hence, transformation to the hyperspherical basis 
could result in violation of the Pauli exclusion principle. 

6. Norm kernel of three a-particles 

Translation-invariant norm kernel for 3a system can be obtained by squaring the norm kernel for three- 
dineutron system: 

/3„ (a,B;A,b) =3! 

For convenience we divided the antisymmetrization operator A of the 3a-system into six parts Ai, each 
corresponding to a certain group of nucleon permutation^^. Six parts /g^ ^A, B; A, B^ of the norm kernel 

Isa ^A, B; A, B^ , which correspond to constituents Ai of the antisymmetrization operator A, are given in 
an explicit form in Appendix |^ 

In this section we are interested in the states of positive parity and zero orbital momentum. Hence, all 
the Pauli-allowed states belong to the (2j/ — 4/i, 2^) representation of the SU(3) group, like in the case of 3d' 
system. So specifying operators Ai we have taken into account the symmetry peculiar to the basis functions 
2/i) (-^' which these operators should act. Below we shall explain this statement in detail. 

Example of the 3d' system illustrates that a permutation of nucleons belonging to different clusters can 
be reduced to the following mathematical operation: 

/ (A, B) ^ const • / (aA + /3B, (3* A + a*B) , (11) 

where a, (3 are some complex numbers. Hence, due to the above-mentioned permutation, some rearrangement 
of the clusters occurs. Then cyclic permutations of such rearranged clusters can be described as 

/ (A, B) const ■ f (e±*^ (aA + /3B), e^'^ {/3*A + a*B)) . 
Obviously, the latter cyclic permutations should give the same result as the initial permutation (jlip : 
^\2l%X) + I3B, l3*A + a*B) + (e^^ (aA + /3B), e^^ (/?* A + a*B)) + 

+^l2.^':X) (^"'^ + A + "*^)) = 3V'(T.-4;'2;.) ("A + /5B, P*A + a*B) . 

In addition, functions 4^2'^^^'^ 2 ) (A-i B) are symmetric with respect to the interchange of vectors A and 
bE] By this reason, the following equality holds true: 

V'(TZTm'2^) + a*B, aA + /3B) = ^'(T.XX) («A + f3B, (3* A + a*B) . 

Keeping in mind all of the preceding, let us determine operators Ai- The first part Ai of the antisymmetriza- 
tion operator A for three a-particles is very similar to the antisymmetrization operator A of three-dineutron 
system: 

^ We denoted an antisymmetrization operator of 3a system by calligraphic A in order to distinguish it from antisymmetrization 
operator A of 3d' system 

^ Such interchange corresponds to the permutation of two boson clusters. 



/ A B A B 



6 

5^-Aexp{(AA) + (BB)} 
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^1 = 6^11 + 6^12 + 3 (A^s + Ms + A3) > 



Ai/(A,B) = /(A,B), ^2/(A,B)=2/(^-^,-|), iJg/CA, B) = 2/ + ^ ^ + ^ 
Operators Af^ are defined as follows: 



AV(A,B)= i°/(A,B 



A^exp(± i^)A,B— >exp(q;i^)B 



Permutations belonging to this group differ from the permutations in three-dineutron system only in one 
respect: the latter refer to fermions, while the former relate to bosons. Considering dineutrons and diprotons 
as bosons, we come to the conclusion that An corresponds to the permutations of a-clusters, A12 can be 
matched with the exchange of dineutrons and diprotons between all a-particles, while Al^ and Af^ are 
associated with the permutations of dineutrons and diprotons between two a-clusters with the third cluster 
being a spectator. 

The second part A2 of the antisymmetrization operator of 3a system 

A2 = 6A21 + 6A22 + 3(^23 + -423 + -423) 
combines permutations involving deuteron clusters: 

i2i/(A, B) = 4/ (-y , -f ) , -422/(A, B) = 8/ |) , A%f{A, B) = 8/ 

^21 includes permutations of deuteron clusters with unit spin projection; A22 is associated with the ex- 
change of nucleons between three deuterons, while the remaining three deuterons (characterized by zero 
spin projection) are spectators; ^23 '^^^a combine permutations resulting in the destruction of one a-cluster 
with the subsequent distribution of two constituent deuterons with unit spin projection between two other 
a-clusters and formation of the cluster from two remaining deuterons with zero spin projections. 
Operator A3 

A3 = 6 (A31 + A32 + A33 + i34 + -^34 + -^34) 

includes permutations of protons, with dineutrons remaining intact: 

A,,m. B) = 4/ 2) ; A./(A,B) = 4/ (i±MA, ; 
^33/,A,B) = 4/('i^A,i±i^BV iU(A,B) = 4/''^*-'' ^'^'^ 



A31 combines simultaneous cyclic permutations of spin-up and spin-down protons; A32 and ^33 contain 
cyclic permutations of spin-down or spin-up protons, accordingly; ^34,^^4 are responsible for the one- 
proton exchange between one cluster and each of the remaining two clusters. 

The next part of antisymmetrization operator A is related to the exchange of spin-down protons between 
two clusters accompanied by the identity or cyclic permutations of dineutrons. 

A4 = 6(^4-j^ -\- A^i -\- A^-^ -\- A^2 -^42 A^2 ^~ "^43 "^43 "^41 ) 

In this case dineutrons and one of diprotons are not rearranged. 

'3A + B A + 3B^ 



^°i/(A,B) = -4/(^: 
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^42/(A, B) = -4/ I ^ , 1 , ^^3/(A, B) = -4/ I , 1 . 

.441, -^42 and Ais are associated with the identity, clockwise or contraclockwise cychc permutations of 
dineutrons, accordingly. 

Operator A5 is characterized by the odd number of pair permutations of nucleons: 

A = 6(3Al + -452 + -4+ + + -^"3 + -4^3 + -^^3)- 

In this case one diproton remains a spectator. 



A,/(A,B) = -8/ (^j.j^ AUiA. B) = -8/ (^Ia + i^e-""=B, -'^e^'/'A 



Ib 

4 



AtJiA, B) ^ -8/ (^\a !^e^W3B, ^^e^^^^A + Jb j . 

A51 combines permutations resulting in the two-nucleon exchange between all clusters; .452, A^2 ^^'^ -^53' -^53 
correspond to the one-nucleon exchange between two clusters accompanied by the exchange of deuterons 
with unit or zero spin projection, accordingly. 

The last part Ae of the antisymmetrization operator looks like 

As - 3(i°i + i+ + Aei + A"e2 + A^^ + A^^) + HA^s + + igs)- 
All constituent parts of operator Ae can be expressed via previously defined permutation operators; 

JO _oJ0 JO —oJO JO _oJ0 

Taking into account that y^2i — 2,4i2 and ,422 = 2^31, we come to the final expression for the antisym- 
metrization operator of the 3a system: 



6 I 

A^Yl T\^i ^ -^11 + ^-^12 + 3^31 + A32 + A33 + 3Ai 



6 

3! 

-^a?3+-4+3+ir3) + ^( 

- Al, + At + A7, + Alo + At + AJ. + At + At + Aj, + At + At + i.„ + At + At + A: 



^(i?3 + i+3 + ir3) + o (-^23 + -4^3 + i^3) + 3(i°4 + i+ + ^34) 



Ml ~r -^41 T- .^41 T- -^42 -T .^42 T .^42 T .^43 T .^43 -r .^43 T .^52 T .^52 T .^52 T .^53 -r .^53 -r .^53. 

Applying operator A to the basis functions ■0|^2T'-^4^^2^i) shall come to the following expansion of 

the partial norm kernel if;^^_4„ 2u) ^'^^ 3a-system: 



*(2"-4m,2m)= Y1 Y1 (2i^,t^-3(7-2/i|2i/,t/-3(7-3q-2/i)3a 



q=0 q=0 
X 



{ J,'^~3q~-2fi Ti/-3g-3g-2/^ . , i/-3g-3g-2p Ti/-3g-2p \ 

\^V^(2i/-4;i,2Ai)^(2i/-4Ai,2Ai) ^ ''^(2i/-4p,2;i) 'f' {2iy~4fj.,2ti} J ' ^-^J 

The above expansion differs from similar expression ([5]) for the threc-dincutron system only in the explicit 
form of the matrix elements of the secular equation: 



{2^, V — 3q ~ 2/i|2i/, 1/ — 3(7 — 3g — 2/i)3Q = 



869-1 V16, 
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m=0 



2m 



iv + 3q + '3q- 2^i)l{iy ~ 3q ~ 3q ~ 2/i)! 



v + 3q — 2^ 



42-1 y + 3q - 2A.)!(i. -3q- 2m)!(1 + (5,,o)(l + ^^,-,) U - 3q - 3g - 2^ 



22^(1 + (-!)■?) 



26«+3Si^i (a, /3; 7; 4) + 22--6g-3g-4,,^ (^a, 7; J 

36«+3,-2 (^^ p. ^. 9) + 32.-69-3g-4p^^^ (^«, 7; i 

(*V3)6^+3^~2Fi(a,/?;7;3) + 3'^-^^(-l)'(zV3)-6^-3SFi («,/3;7; J) }) 



where indices of hypergeometric function are defined as follows: 

a = —V + 3q + 2/i, /3 = —V + 35 + 3(7 + 2/i, 7 = 6g + 3(2" + 1. 

Eigenfunctions ^'(21^-4;^. 2^1) of the antisymnietrization operator A for the 3a-system are given in terms of 
orthogonal polynomials of a discrete variable, as in the case of the 3d'-system (see Eq. ([7])). The number 
of different states possessing SU(3) symmetry (2v — 4/i,2/i), i.e., the degree of the SU(3)-degeneracy, is 
the same for three-boson system, 3d'-system and 3Q!-system and equal to \"~^^ \ + 1- However, some of 
these eigenfunctions correspond to zero eigenvalues and thus they are forbidden by the Pauli principle. The 
number of the latter functions depends on the system considered and increases with increasing the number of 
nucleons in each cluster. For example, in the three-boson system all the states belonging to the (2i^ — 4/j,, 2/i) 
representation have been allowed by the Pauli principle, while in the three-dineutron system at a given 
number of oscillator quanta v there is a Pauli-forbidden state ^(2jy.o)fc=o- the 3Q!-system at > 3 there 
are two Pauli-forbidden states belonging to the (2i/, 0) SU(3)-representation and one Pauli-forbidden state 
characterized by the (2v — 4, 2) SU(3)-indices. All the Pauli-allowed states of the 3a-system have distinct 

Table 1 

Eigenvalues A(2^_4fi,2n)fc belonging to the first three families of the Pauli-allowed states of the 3o-system 





2fc = 4 




2fc = 6 




2fc = 8 


V 


fj. = 2 fi = 1 fi = 


M = 3 p 


= 2 11 = 1 /i = 




= 3 fi = 2 fi = 1 fj, = 


4 


0.5933 










5 


0.6748 0.4894 










6 


0.7161 0.5647 0.6497 


0.8714 








7 


0.7366 0.6333 0.6269 


0.9081 


8657 






8 


0.7450 0.6859 0.6457 


0.9246 


8852 0.9073 


0.9640 




9 


0.7482 0.7196 0.6773 


0.9324 


9042 0.8959 0.9478 


0.9754 


9666 


10 


0.7493 0.7372 0.7071 


0.9356 


9191 0.9000 0.9193 


0.9805 


9712 0.9770 


11 


0.7497 0.7449 0.7283 


0.9368 


9287 0.9106 0.9084 


0.9828 


9759 0.9738 0.9865 


12 


0.7499 0.7480 0.7404 


0.9372 


9338 0.9212 0.9093 


0.9838 


9796 0.9747 0.9792 0.9906 


13 


0.7500 0.7492 0.7461 


0.9374 


9360 0.9290 0.9159 


0.9842 


9821 0.9774 0.9763 0.9837 


14 


0.7500 0.7497 0.7485 


0.9375 


9369 0.9336 0.9236 


0.9843 


9834 0.9801 0.9764 0.9792 


15 


0.7500 0.7499 0.7494 


0.9375 


9373 0.9358 0.9298 


0.9843 


9840 0.9821 0.9783 0.9773 


16 


0.7500 0.7499 0.7498 


0.9375 


9374 0.9368 0.9337 


0.9844 


9842 0.9833 0.9804 0.9775 


17 


0.7500 0.7500 0.7499 


0.9375 


9375 0.9372 0.9358 


0.9844 


9843 0.9839 0.9822 0.9790 


18 


0.7500 0.7500 0.7500 


0.9375 


9375 0.9374 0.9368 


0.9844 


9844 0.9842 0.9833 0.9808 


19 


0.7500 0.7500 0.7500 


0.9375 


9375 0.9375 0.9372 


0.9844 


9844 0.9843 0.9839 0.9823 


20 


0.7500 0.7500 0.7500 


0.9375 


9375 0.9375 0.9374 


0.9844 


9844 0.9843 0.9842 0.9833 
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eigenvalues, which tend to the eigenvalues of a 2a subsystem with increasing the number of oscillator quanta 
u (see Table [J): 



lim A 



1^ — 2/^- 



(2u-4tJ..2fi)k 



\^Bo=2a 



fc-1 



k,0- 



Eigenfunctions ^(2v-4iJ,,2fj,)k ^1^° become more complicated with increasing the number of nucleons. In 
the case of Sa-system each eigenfunction of the antisymmetrization operator is a superposition of func- 



tions i^'(2,%^.2^,) 



with all nonvanishing expansion coefficients D 



l^ — 3q—2fj, 
(2i/-4/i,2M)fe- 



However, as the eigenvalues 



A 



(2i/-4^,2Ai)fe 



of the 3a norm kernel approach limit values A2^°, corresponding eigenvectors (2iy-4fj..2ti)k take 



simple analytical form 5'?2i/_ 



(2i/-4^,2M)fc 



(see Eq. ([9])), as can be seen from Tabled Consequently, in the limit 

ly ^ k the expansion coefficients D^(2u~i^2ti)k '^^'^ identified with the Kravchuk polynomials of a discrete 
variable: 



3g— 2/i 
^(2v-4n,2ii)k 



1 + 5, 



9,0 



^2^ (^-3gr-2/i), 



where the explicit form of the Kravchuk polynomials /Cjj, iy — iq — 2/i) is given by Eq. ([11 

Therefore, index fc, which makes sense of the number of oscillator quanta accounted for by the ®Be sub- 
system, serves as an additional quantum number for the SU(3)-degenerated states. It should be noted also 
that due to the incoincidence of eigenvalues ^(^-^^1^ 2^1)^^ eigenfunctions ^'(2i/-4^,2/j)fc are uniquely deter- 
mined. Any unitary transformation applied to the latter SU(3)-basis functions would disrupt the diagonal 
form of the norm kernel of the 3a system and, hence, is inappropriate in this case. 



Noteworthy also is a remarkably simple form of the eigenfunctions [2v- 



y^(2i/-4M,2M)fc=j, - V-'(2,y-4Ai,2M) ' 



— 2/i 



= 0. 



ifi,2^i)k=^ 



in the case — 2/i < 3 : 



It is precisely this fact that explains why the overlap integrals J ^ {2v-i^^,2tJ,)k=^l^\2l'-'^^^ 2^i)k- '^'^^ ^''^ equal 
to unity, as long as z/ — 2// < 3. 

In turn, the eigenvalues ^{2v-i^j..2ti.)k= are given by the following analytical expression, provided that 
J/ - 2/t < 3 : 



A 



(2iy-4/i,2^)fc=,, 



6 

4^ 



1 



2 

4^ 



8 



18 
42I7 



2Fi{-u + 2/i, -V + 2/i; 1; 4) + 4^-2^2^^^ + 2/i, ~v + 2/i; 1: 
2Fi{-v + 2/i, -^v + 2/i; 1; 9) + O'^-^^ai^i ( + 2/4, + 2/t; 1 



2Fi{~v + 2/i, ~v + 2/i; 1; 3) + r-^^'2Fl -v + 2/i, -v + 2/i; 1 



32^ 



i82M 

3 

22^+1 



From Table[l]and Tablc[2]we notice that the Pauli-allowed basis states 5'(2i/-4/i. 2/^)1- can be arranged into 
the branches and families, with all the states of a particular branch having common SU(3)-symmetry index 
\i and overlapping generously with corresponding asymptotic function ([HI), but differing in value of the first 
index A of the SU(3)-symmetry. The eigenvalues belonging to a given branch tend to the same limit value 
^2k^ with the number of quanta increasing. The branches which share limit eigenvalues are combined in 
the family of the eigenstates, which thus is completely determined by the degree k of the corresponding 
Kravchuk polynomial. 

It should be noted that the states listed in Table [T] and Table [2] exhaust all possible basis functions allowed 
by the Pauli principle, as long as < 9. As increases, new families of states characterized by 2fc > 8 also 
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Table 2 

Overlap integrals / ^(2i^-4/j.2;t)fc '^'^2i'-4p 2p) , '^'^s versus the number of quanta v for the first three families of the Pauli- 
allowed states of the So-system 
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1 
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0000 
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appear. However, eigenstates and eigenvalues belonging to such families are governed by the same laws as 
those tabulated in Tables [HO The state ^(o,4)j._2 characterized by the minimum number of quanta u = A 
and belonging to the fi — 2 branch of the k ~ kniin = 2 family of the Pauli-allowed states is identical to 
the state of the leading SU(3)-representation of the oscillator translation-invariant shell model in the Elliott 
scheme. Other states belonging to the /i = 2 branch provide the main corrections to the states of the Elliott 
model resulting from the a-clustering of the ^^C nucleus and its decay through the ^^C— >®Be-|-a channel. 
Branches of the states belonging to the families characterized hy k > k^i^ reproduce the excitations of the 
^Be accompanied by the increasing of the number of oscillator quanta accounted for by this subsystem. 

Existence of different limit eigenvalues A2^° reflects the possibility for two a-clusters to be close to each 
other and far apart from the third a-cluster. Owing to this the unitary transformation from the SU(3)- 
basis to the other one is rendered possible only within a particular family of the Pauli-allowed states. 
Transformation to the three-cluster hyperspherical harmonics, contrastingly, involves the states belonging to 
different families. The hypermomentum can not serve as quantum number of Pauli-allowed basis functions, 
because the Pauli exclusion principle mixes the states with different values of hypermomentum. Hence, 
restriction for the maximum value of hypermomentum K included in calculation K < K„iax leads to 
spoiling the Pauli-allowed basis functions corresponding to the number of oscillator quanta v > K„iax , and 
this effect is enhanced with increasing v. The angular-momentum-coupled basis can be used instead in the 
region v ^ k, where all eigenvalues ^(2iy-A^i,2ii)k belonging to the fcth family coincide with their limit values 
•^2fc°- Basis functions of the latter basis are labeled by the number of quanta, and the angular momenta 
of the ^Be subsystem coinciding with the angular momentum of the relative motion of this subsystem and 
a remaining g-particlj^. The transformation to this basis is performed within a particular family of the 
Pauli-allowed states and, therefore, does not result in violation of the Pauli principle. At the same time, the 
angular-momentum-coupled basis allows one to set the asymptotic boundary conditions for the expansion 

* Recall that only the states with total angular momentum L = are considered in this paper. 
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coefficients of a three-cluster wave function in the continuum. This point has been discussed in detail in Ref. 
[13] with the •^H+n + n system. 

From the above discussion it appears that asymptotically each family of Pauli-allowed states of Sa-system 
corresponds to a certain binary decay channel of the latter system into a 2a subsystem occurring in a 
ground or an excited harmonic-oscillator state and a remaining a-particle. Hence, the excited states of the 
^^C nucleus should decay via subsequent stage ^^C^^Be+a a + a + a rather than in a "democratic" way 
-'^^C^ a + a + a usually associated with hyperspherical harmonics. The latter basis can, probably, provide an 
adequate description of the ground state of the ^^C nucleus, but it is definitely not appropriate for studying 
the continuum states. 

Finally, let us compare the behavior of the eigenvalues and eigenfunctions of the Sa system with those 
of the ^H='^H-|-n -|- n system discussed in detail in Ref. [l^. First of all, at a given number of quanta i/ 
the number of the Pauli-allowed states in the 3a system is much less than that of the system. In the 
3a system only the families of the states corresponding to even values of quantum number k are consistent 
with the requirements of the Pauli exclusion principle, while the model space of the Pauli-allowed states of 
'^H+n + n system includes the families characterized by both odd and even values of k. The number of the 
Pauli-forbidden states in the 3a-system is also greater than in the nucleus. In the system only one 
branch falls into the category of the Pauli-forbidden states. It is characterized by fc = and corresponds to 
the occupation of an s-state in the '^H subsystem by one of the valence neutrons. The 3a system have three 
branches of the Pauli-allowed states: a branch being a member of fc = family and two branches entering 
in the family 2fc = 2. These two families correspond to two forbidden states of the 2a subsystem. Moreover, 
whereas in the case of the nucleus all branches belonging to the same family appear at the same number of 
quanta, the branches combining into a particular family of the 3a system start in one-quantum increment. A 
new family of the Pauli-allowed states in the system arises with increasing the number of quanta by one, 
while families of such states in the 3a system emerge in three-quantum interval. This results in decreasing 
the model space of the Pauli-allowed functions in the internal region of small distance between clusters in 
the 3a system as compared with the system. 

All the differences discussed above are due to a high symmetry of the 3a-system, which comprises of 
three identical clusters, while the ^H+n + n system have only two identical constituents. In addition, in 
the 3a-system each family begins with the branch having maximum possible SU(3)-symmetry index fi — k. 
This fact gives grounds to expect that the branches characterized by = fc should be more favorable than 
the branches with /i = 0, in contrast to the case of the nucleus. Of course, both for the case of the 3a 
system and the "^H+n + n system the families of the allowed states associated with the low-order Kravchuk 
polynomials should dominate. 

The eigenvalues of the antisymmetrization operator of the 3a system tend to unity from below resulting 
in the repulsive effective interaction between a-clusters due to the change in the kinetic energy of relative 
motion under the effect of the Pauli exclusion principle (see Refs. |13|19j for details). Hence taking into 
account all the exchange effects in 3a-system may allows one to overcome the problem of overbinding of the 
^^C nucleus without introducing an additional three-cluster repulsion. 

7. Conclusion 

Within the microscopic model based on the algebraic version of the resonating group method the role 
of the Pauli principle in the formation of continuum wave function of nuclear systems composed of three 
identical s-clusters has been investigated. Our principal concern has been with the study of the exchange 
effects contained in the genuine three-cluster norm kernel, i.e., taking into account the eigenvalues of the 
Pauli-allowed states. 

Norm kernels for the system of three identical fermions and the system of three identical bosons have been 
constructed in the Fock-Bargmann space. It has been shown that the former serves as the main building 
block for the norm kernels of nuclear systems composed of three identical clusters, while the latter generates 
the basis functions needed for the description of nuclear systems composed of three identical boson clusters 
(such as n +'^ n +'^ n and ^^C=a + a + a). Simple analytical method of constructing the norm kernel 
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for 3d' system and 3a system in the Fock-Bargmann space has been suggested and realized. 

The Fock-Bargmann image of the antisymmetrization operator has been found for three-fermion, three- 
boson, three-dineutron and 3a systems. Classification of different constituents of the antisymmetrization 
operator by belonging to a certain group of permutations is given. Careful analysis of the structure of the 
eigenfunctions and behavior of the eigenvalues of the antisymmetrizer has been performed for all the above- 
mentioned three-cluster systems. It has been demonstrated that the Pauli-allowed states of three-cluster 
systems composed of identical boson clusters can be represented as a superposition of the eigenfunctions of 
the three-boson norm kernel, with the expansion coefficients being the orthogonal polynomials of a discrete 
variable. 

The Pauli-allowed basis functions for the 3a and the 3d' systems are given in an explicit form and the 
asymptotic behavior of these functions is established. Eigenvalues of the three-cluster system are shown to 
tend to the eigenvalues of a two-cluster subsystem with increasing the number of oscillator quanta v. At 
the same time, the corresponding eigenvectors take simple analytical form, while the expansion coefficients 
become the Kravchuk polynomials as the number of oscillator quanta increases. We suggest a way of resolving 
the problem of the SU(3) degeneracy of the Pauli-allowed states. A degree of the Kravchuk polynomial serves 
as an additional quantum number to label the states belonging to the same SU(3) irreducible representations. 

Complete classification of the eigenfunctions and the eigenvalues of the ^^C norm kernel by the ^Be= 
a + a eigenvalues has been given for the first time. We have demonstrated that for the 3a system such 
classification is unique in that it is consistent with the requirements of the Pauli exclusion principle both 
in the region of small intercluster distances and in the asymptotic region where the scattering matrix 
elements are produced. Due to incoincidence of eigenvalues of the antisymmetrization operator of the 3a 
system the corresponding eigenfunctions are uniquely determined. Any unitary transformation applied to 
the latter SU(3)-basis functions would disrupt the diagonal form of the norm kernel of 3a system and, 
hence, is inappropriate in this case. Contrastingly, the eigenfunctions of the three-neutron or the three- 
boson systems are determined only to the unitary transformation, because all the Pauli-allowed states have 
unit eigenvalues. Hence, in this case the SU(3)-classification is only one of possible variants. The 3d' system 
is a borderline case, because the eigenvalues of the norm kernel are equal and tend to unity for all the 
Pauli-allowed states. Consequently, any transformation, which involves only the Pauli-allowed states and 
remains the Pauli- for bidden states intact, does not disrupt the diagonal form of the norm kernel. 

Summarizing, the Pauli-allowed states of the 3a system can be arranged into the branches and the families. 
The cigcinvalues belonging to a given branch tend to the same limit value Ajf of the 2a subsystem with 
the number of quanta increasing. The branches which share limit eigenvalues are combined in the family of 
the eigenstates, which asymptotically corresponds to a certain binary decay channel of the 3a system into 
a 2a subsystem occurring in a particular harmonic-oscillator state and a remaining a-particle. Hence, the 
excited states of the ^^C nucleus should decay via the subsequent stage ^^C^*Be-|-a a -H a -|- a rather 
than in a "democratic" way ^^C^ a -I- a -I- a usually associated with hyperspherical harmonics. 

Propc;r truncation of the model space of the Pauli-allowed basis functions consists in neglecting families of 
the states with high values of additional quantum number k and using the angular-momentum-coupled basis 
in the region v ^ k for setting the boudary conditions for the three-cluster wave function in the continuum. 

Spectrum of the '^^C norm kernel has also been compared to that of the ^H system. The model space of 
the Pauli-allowed functions in the internal region of small distance between clusters in 3a system has been 
shown to decrease essentially as compared with the ^H system due to a high symmetry of 3a-system, which 
comprises of three identical clusters. 

Appendix A. Partial norm kernels of the 3a system 

The norm kernel of the 3a system can be divided into six partial norm kernels, each corresponding to a 
certain group of nucleon permutations: 

6 6 ^ 

J3a (a, B; a, b) = ^ -A^ exp | (AA) + (BB)} = ^ -/« (a, B; A, b) , 

i=l i=l 
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where 

4a = -^{Xl + yi + Zi + X2 + y2 + Z2f, 1^2 = ^(^1 + Vl + Zl + X2 + y2 + Z2f 

/Q\ 1 

Ha = ^2(^1 +yi+ Zl+X2+y2 + Z2){XI + j/l + ^1 + X2 + ^2 + Z2) 
4a = ^■^^(^l +yi + Zl+X2+y2 + Z2){X3 + ys + Z3){X4 + 1/4 + Z4) 
4a = "■^'^(^1 +yi+Zl+X2+y2 + Z2){X3 + ^3 + Z3){x4 + ^4 + Z4) 

4a = ^4(a;3 + 1/3 + ^;3)^(a;4 + y4 + Z4f. 

We used the following notations 

f(AA) (BB)1 f(BA) (AB)1 

a:i=exp|^ + ^|; -2 = exp|^ + ^|; 

f ,2^ (AA) _,2^ (BB) 1 f (BA) _.2^ (AB) 1 

yi = exp <^ 3 ^—1 + e ^ 3 ^ , = exp <^ 3 ^—1 + g * 3 ^—^ ^ 

f _,2^(AA) i2.(BB)l f „,2,(BA) .2,(AB)1 

= exp <^ e ^3 ^--2 + 3 ^ Z2 = exp <^ e * 3 ^—1 + 3 \ 

_ ^ / i2^(AA) _.2^(BB)1 _ ^ / i2^(BA) _i2^(AB) 
1/1 = 2 exp < -e' 3 ^ ^ ' - e ' 3 -L__^ ^ 1/2 = 2 exp < -e' 3 ^ ^ ^ - e ' 3 -L__^ 

_ ^ / _i2^(AA) i2^(BB)l _ ^ / _i2^(BA) ,2^(AB) 
2:1 = 2 exp < -e * 3 ^ ^ ^ - e' 3 k ^ Z2 = 2 exp < -e ' 3 ^ ^ ' - e' 3 -L__^ 

f i2,(AA) _,2,(BB)1 f i2,(BA) _,2,(AB)1 

y3 = exp I 3 + e ' 3 ^ , t/4 = exp ^ 3 + e ' 3 V 

f _i2^(AA) i2^(BB)l f _i2^(BA) i2^(AB)l 

23 = exp < e '3 ^ ^ + e* 3 ^—-L K 24 = exp < e * 3 ^—-J- + 3 ^ V 
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